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Abstract. - We study random walks on a family of treelike regular fractals with a trap fixed on a 
central node. We obtain all the eigenvalues and their corresponding multiplicities for the associated 
stochastic master equation, with the eigenvalues being provided through an explicit recursive 
relation. We also evaluate the smallest eigenvalue and show that its reciprocal is approximately 
equal to the mean trapping time. We expect that our technique can also be adapted to other 
regular fractals with treelike structures. 



Introduction. — As a valuable tool for various scien- 
, tific disciplines, fractals introduced by Mandelbrot [T] have 
' received considerable attention in many fields in the past 
, few decades [2]-[5] . The so-called regular fractals constitute 
' an integral family of fractals, including the Sierpinski tri- 
, angle [6] , hierarchical lattice fractal [7] , Vicsek fractals [8] , 
' T fractal ^^U\ , Peano basin fractal [TT] , to name a few. 
, These simple fractal structures have been a focus of ex- 
' tensive research [4] . A great advantage of regular fractals 
, is that many problems about them can be treated analyt- 
' ically, presenting a first step in understanding the struc- 
, tural and dynamical proprieties of general fractals. 
] A fundamental problem in the study of fractals is to find 
I the relations between structural properties of fractals and 
' their dynamics, which remains one of the great challenges 
1 in the realm of nonlinear science. As a paradigmatic dy- 
namical process, random walks on fractals play a crucial 
role in a wide range of applications [2111], therefore have 
been intensely investigated [T2U22j . These investigations 
have uncovered some nontrivial effects of topological struc- 
tures of different fractals on the behavior of random walks. 
Among various random- walk issues [23l[2l|, the trapping 
problem on fractals are becoming increasingly important 
since many phenomena built on fractals can be described 
by random walks [251 [26] . Particularly, the mean first- 
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passage time (MFPT) to the trap node in a network has 
been widely studied [?7H5i] due to its role as an indicator 
of the transport efficiency. 

In view of the importance of MFPT [35] , a huge amount 
of effort has been devoted to developing techniques to 
determine MFPT to a given node, and numerous meth- 
ods have been proposed, including the generating func- 
tion [27], renormalization group [28l|3T], back- forward 
equation [35], and so on. Most of existing methods are 
presented for particular networks and thus are applicable 
only to some specific structures. A universal technique 
for general networks is still unavailable, with the excep- 
tion of the spectral approach, through which the MFPT 
from a node to another can be expressed in terms of the 
eigenvalues and their corresponding eigenvectors of a re- 
lated transition matrix |36j . For the trapping problem, 
the average trapping time (ATT), i.e., the MFPT to the 
trap, can be approximated by the inverse of the smallest 
eigenvalue of an associated matrix [37]. Irrespectively of 
the generality of the spectral approach, the computational 
complexity for determining the spectra of generic networks 
is very high [35], which makes the task very difficult and 
even impossible for large networks. In this context, it is 
of interest to seek new techniques for evaluating the spec- 
trum of a relevant matrix for some particular networks. 

In this paper, we study a random walk problem on a 
family of treelike regular fractals [39] with a perfect trap 
fixed at a central node. The fractals under investigation 
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Fig. 1: (Color online) Construction of the fractals. The next 
generation is obtained through replacing each edge by the clus- 
ters on the right-hand side of the arrow. 



include both the T fractal [HIII^ and the Peano basin 
fractal [TT] as their particular limiting cases. Using the 
real-space decimation procedure, we derive an explicit re- 
cursive relation for the eigenvalues of a matrix associated 
with the trapping issue. We then continue to determine 
the full spectrum of eigenvalues and their multiplicities 
of the relevant matrix. In the end, we derive a recursion 
relation for the smallest eigenvalues at two consecutive 
iterations of the fractals, based on which we obtain an ap- 
proximate value of the final smallest eigenvalue and then 
further evaluate the ATT. 

Constructions and properties of regular fractals. 

— The fractals under consideration are constructed in 
the following iterative way l39|. Let Fg {g > 0) denote the 
fractals after g iterations. Initially {g = 0), Fq consists of 
an edge linking two nodes. In each subsequent iteration 
g > 1, Fg is obtained from -Fg-i by performing the fol- 
lowing operations on each existing edge as illustrated by 
Fig- HI replace the edge by a path two-links long, with the 
two end points of the path being the same end points of 
the original edge; then, create m (an natural number) new 
nodes and attach them individually to the middle node of 
the path. Figure [2] describes the first three steps of the 
iterative construction process for a particular fractal cor- 
responding to TO = 1. The family of fractals, controlled 
by parameter m, includes several well-known fractals as 
its limiting cases: When m = 1, it reduces to the T frac- 
tal [SlIlO]; when m = 2, it recovers the Peano basin frac- 
tal [TT]- 

The studied fractals have an obvious self-similar struc- 
ture, which can be seen from their alternative construction 
approach. If we call the central node (node 3 in Fig. [21 
for example) the innermost node and those nodes farthest 
from the central node the outermost nodes, then the frac- 
tals can also be built as follows (see Fig. [3]): Given the 
generation g, we can obtain -Fg+i by merging TO-f 2 copies 
of Fg with the to -I- 2 outermost nodes in separate replicas 
being coalesced into a single new node, i.e., the innermost 
node in Fg^i. 

It is easy to derive that the total numbers of edges and 
nodes in Fg are Eg = [m + 2)» and Ng = {■m + 2y + 1, 
respectively. Some other properties of the fractals can 
also be determined [S^. Their fractal dimension is df = 
ln(TO+2)/ In 2, while their random- walk dimension is = 




Fig. 2: (Color online) Iterative process for a special fractal 
corresponding to m = 1. 

ln[2(TO + 2)]/ In 2. Thus, their spectral dimension is = 
2di/d^ ^ 21n(TO + 2)/ln[2(TO + 2)]. 

Random walks with a trap positioned on the in- 
nermost node. — In this section, we study the discrete- 
time unbiased random walks taking place on the fractals 
Fg. In the process of the classical random walks, at each 
time step, the particle (walker) starts from its current lo- 
cation and moves to any of its nearest-neighboring nodes 
with the same probability equal to the reciprocal of the 
degree of the current node. 

Now, consider the probability distribution function 
Pj [t) , which is the probability for the walker being on node 
J at time t. If the starting point for the walker is s, then 
the stochastic master equation governing the evolution of 

Pj{t) is BOlllI] 

^--E^.P.w (1) 
i— 1 

with the initial condition dpj{t = 0) = Sj^s- 

In Eq. ([T|), if i is not equal to j, is the jumping 
probability from node i to node j, but Pij = — 1 other- 
wise (i = j). Note that we can consider as an entry 
of matrix Pg, called the transition matrix, which suggests 
that the powerful methods of the spectral theory for ma- 
trices [12] can be applied. 

Eigenvalues of the stochastic master equation. We fo- 
cus on a trapping problem which is a particular case of 
random walks on Fg with a trap fixed on the central node 
of the fractals. For this case, the transition matrix is a 
submatrix of Pg, which is obtained from Pg by deleting 
the row and column corresponding to the trap node. In 
the case without any confusion, we also use Pg to repre- 
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Fig. 3: (Color online) Second construction method of the frac- 
tals highlighting their self-similarity. We can obtain -Fg+i by 
joining m + 2 copies of Fg, represented by Fg^\ Fg^\ . . ., 



and 

an outermost node. 



X is the innermost node, while K is 



sent the transition matrix of the trapping problem. Below 
we will use the decimation method [43H45] to determine 
the complete spectrum of eigenvalues and their degenera- 
cies, and demonstrate that the reciprocal of the small- 
est eigenvalue has the same scaling as that of the ATT. 
Notice that previous works have studied analytically the 
Laplacian spectra for different treelike graphs, e.g., den- 
drimers [46], Vicsek fractals [47l|48], and related struc- 
tures, say Husimi-cacti [33] ■ Here we concentrate on the 
normalized Laplacian matrix, since for an irregular graph, 
its Laplacian matrix and normalized Laplacian matrix be- 
have quite differently, and the latter seems to be more 
natural [50] . 

The decimation method has an advantage that one can 
solve the eigenvalue equation based on a similar problem 
from the previous iteration. Let a denote the set of nodes 
belonging to Fg, and /3 the set of nodes generated at iter- 
ation 5 + 1. Then, one can write the eigenvalue problem 
for the transition matrix Pg+i in the following block form: 



P. 



P, 



Po^.P 



Up 



= A 



3+1 



Up 



(2) 



where Pa, a is the identity matrix since there are no tran- 
sitions between the original nodes, Ppjj is block diagonal 
with each block being the same (m -I- 1) x (m -1-1) matrix 
of the form 



B = 



^ m+2 m+2 



-10 

and \g+i is an eigenvalue of Pg+i- 



(3) 



Equation ^ can be rewritten as two equations: 

Pa,aUa + Pa,l3Ufj = Ag+iWa, (4) 
P0,aUa + Pl3,pUp = Xg+lUp. (5) 

From Eq. ([S]), one obtains 

Uf} = (Ag+i - Pp,l3)~^Pp,aUa , (6) 

provided that the concerned matrix is invertible. Inserting 
Eq. © into Eq. @i leads to 

[Pa,a + Pa.pi^'g+l ^ Ppj)'^ Pp.alUa = Xg+lUa- (7) 

Thus, the problem of determining the eigenvalue Xg+i 
for matrix Pg+i of order X is reduced 

to computing the eigenvalue problem of a matrix with a 
smaller order of (m + 2)^ x (m + 2)^. 

Let Qg = [Pa,a + Pa,l3{Xg+l - -P/3,/3 ) " ^ -P/3,a] ■ Ncxt, We 

will show that matrix Qg can be related to matrix Pg, and 
thus the eigenvalue A^+i for matrix Pg+i can be expressed 
in terms of eigenvalue Xg for matrix Pg. To this end, we 
define 

R.g = {x + y)Ig-yPg, (8) 
where Ig is the identity matrix of the same order as Pg, 



1 + 



Xg+l - 1 



and 



(m + 2)(Ag+i -2)Ag+i +2' 
A„ 



y 



1 



(m + 2)(Ag+i -2)Ag+i +2' 



(9) 



(10) 



Equation ^ indicates that the eigenvalues of Rg are re- 
lated to those of Pg. In the Appendix, we prove that 
Qg = Rg, thus the eigenvalues of Pg+i can be obtained 
through determining the spectrum of Pg . 
Substituting Rg for Qg, Eq. ([7]) becomes 



which means 



that is. 



= {x + y)IgUa - yPgUa 

= Xg+lUa, 



{x + y- Xg+i)Ua = yPgUa, 



y 



y 



(11) 



(12) 



(13) 



If Ag denotes the eigenvalue of Pg corresponding to the 
eigenvector Ua, Eq. ([13]) implies 



A„ 



y 



X 



9+1 



(14) 



Plugging Eqs. ([9]) and ([T0| in Eq. ([M]) . one can see that 
the eigenvalue Ag+i for transition matrix Pg+i is related 
to that of Pg via 



'(TO + 2)(Ag+l -2)Ag + l, 



(15) 
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which can be recast as 

\2 



(m + 2) {\g+iY ~ 2(m + 2)Ag+i + A,, = 0. (16) 
Solving Eq. gives its two roots 



A 



9+1 



= l± \ l- 



A„ 



2 ' 



(17) 



which provide two eigenvalues of Pg+i- 

By calculating the eigenvalues numerically, one can find 
some important properties of the spectra. First, all eigen- 
values appearing at a given generation gi always exist at 
the next generation + 1. Second, all new eigenvalues at 
generation gi + 1 are just those generated via Eq. ([T7|) by 
substituting Ag with Ag^ that are newly added to gener- 
ation Qi. That is, all eigenvalues can be obtained by the 
recurrence relation Eq. ([T7]). Then, the next step is to de- 
termine the degeneracies of the eigenvalues, based on the 
two characteristics. 

Multiplicities of eigenvalues. Let Ng{X) denote the 
multiplicity of eigenvalue A of the transition matrix Pg. 
Since each eigenvalue in the fractals can be derived from 
eigenvalue 1, we first determine the number of eigenvalue 
1 of Pg. For this purpose, let r{M) be the rank of matrix 
M . Then, the number of eigenvalue of Pg is 



iVg(A = 1) = {m + 2Y 



riPg 



1 X /„ 



(18) 



For the two limiting cases oi g = 1 and g = 2, we can 
easily obtain iVi(A = 1) and A^2(A = 1). When g = 1, 
Pg is the identity matrix of order (171 + 2) x (m + 2), so 
Ni{X = 1) = m + 2. When g = 2, Pg - Ig is block 
diagonal, with each of its m + 2 blocks being the same 
(m + 2) X (m + 2) matrix in the form of 





-1 
-1 



1 1 

m+2 m+2 



1 

' m+2 










(19) 



Obviously, the rank of each block is 2. Thus, we have 
7V2(A = 1) = (m + 2)2 - 2(?7i + 2) = m(m + 2). 

In order to determine Ng{X = 1) for the case oi g > 2, 
we define Bg as the transition matrix of Fg with the trap 
located at an outmost node, e.g., node K in Fig. [3l and 



define Ag = Bg 



We next show that it is possible to determine exactly 
the multiplicity of eigenvalue 1 for the fractals at every 
iteration g > 2. 



Note that for g > 2, Pg^ 




+1 has the following 


form: 




" Ag 
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(20) 
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where Ag (g > 2) satisfies 



A,= 



Ag-l 








■ -Vi 





Ag-l 





■ -V2 








Ag-l ■ 


■ -V3 


-v7 






■ Ag-l 



(21) 



in which each Vi (i = 1, 2, • ■ • , m 4- 2) is a matrix of order 
(m + 2)f X (m + 2)3, which has (m 4- 2)^^ — 1 elements 
and only one nonzero element l/(m -I- 2) describing the 

. (i) 

transition probability from one node in a replica of Fg 
to the innermost node being joined together to form Pg+i 
(see Fig. [3]); and the superscript T of a matrix represents 
its transpose. 

By construction, for any node u in Fg{g > 2) that is di- 
rectly connected to the center node, it must have a neigh- 
bor h with degree 1. Thus, in matrix Ag, there is only one 
nonzero element for row h and for column h, respectively, 
namely, {h,u) ~ —1 and {u,h) = — l/(m-f 2). Hence, 
we can use elementary operations of a matrix to eliminate 
all nonzero element —l/{m + 2) on the last row and the 
last column of Ag. Then, wc have r{Ag) = r{Pg — Ig), 
which imphes that r{Pg+i — Ig+i) = (m + 2)r[Pg — Ig). 
This, together with r{P2 — h) = 2(m + 2), leads to 
r{Pg — Ig) = 2{m + 2)3~^. Combining the above-obtained 
results, the degeneracy of eigenvalue 1 is 

Ng{X = 1) 



m + 2, 

m{ni + 2)3-^, 



.9 = 1, 

g > 2. 



(22) 



We proceed to determine the multiplicities of other 
eigenvalues. For the convenience of description, on the 
basis of Eq. ^5)) , we define a set 



{1}, 



{x\{m + 2){2x-x^) E^g-i}, g > 1. 



(23) 



It is easy to verify that the cordiality (number of elements 
of a set) of fig is jfigj = 2^"^. On the other hand, from 
the above discussion, we can conclude that all the eigen- 
values of Pg belong to the set {ni,n2,--- ,^g}. More- 
over, extensive numerical computation on the eigenval- 
ues of the studied fractals shows that the multiplicities 
of eigenvalues for the transmission matrix of the fractals 
satisfy A^g(A S ili) = 7Vg_i(A G fli-i). Using this prop- 
erty, we can obtain the degeneracy of each eigenvalue A in 
ri, as 



Ng{X e n,) 



m + 2, 
m(m 4- 2)3' 



« = 5, 
i < g. 



(24) 



Taking into account the degeneracies so obtained, the 
number of eigenvalues of Pg is found to be 



J2Ng{Xen^)\n^ 



1=1 



S-1 

E 



m(m + 2)9-''2*-i + (m + 2)2^ 
{m + 2)3, (2 
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which indicates that we have found all the eigenvalues of 

Smallest eigenvalue and average trapping time. The 
probability distribution function p(t) considered above is 
very important since it describes the evolution of the sys- 
tem. Many quantities of interest are encoded in this un- 
derlying probability distribution function. For example, 
in the limit of large system sizes, the ATT denoted by Tg 
for the trapping problem on Fg is asymptotically equal to 
the inverse of the smallest nonzero eigenvalue of Pg, de- 
noted by Zg, that is, Tg ~ l/^g- Here, Tg is defined as the 
average of trapping time Tg{v) of node v over all nodes ex- 
cept the innermost trap node. Conceptually, Tg{v) is the 
expected time for a walker starting from node v to reach 
the trap for the first time. 

According to the recurrence relation given by Eq. p7)) 
that dominates the eigenvalues of the transmission ma- 
trix for the fractals at two consecutive iterations, we can 
evaluate the smallest nonzero eigenvalue Zg of Pg. From 
Eq. (jl7p . wc obtain the following obvious recursive rela- 
tion: 

' (26) 



Zg+l = 1-^1- 

Using Taylor's formula, we have 



Zg+l W 1 



1 



1 z„ 



2 m + 2 



2(m + 2) 



(27) 



Considering the initial condition zi = 1, Eq. (|27p can be 
solved by induction to yield 



1 



^ ~ [2(m-h2)]9-i 
Then, for large networks. 



1 



2S-i(m-t- 2)3-1 



(28) 



(29) 



We continue to represent Tg as a function of the net- 
work size Ng, in order to derive the scaling between these 
two quantities. From Ng = {m + 2)^ -I- 1, we have 
g = log„^2(-^s ~ Hence, for large systems, Tg can 
be expressed in terms of Ng as 



{T)g ~ (Arg)l+ln2/ln(m+2) ^ (Ngf/ds 



(30) 



Thus, for the entire family of the studied fractals, in the 
limit of large system sizes, the ATT grows as a power-law 
function of the network size Ng. We note that the obtained 
exponent 2/ds perfectly agree with the previously obtained 
result [HlaillMl. 



Conclusions. — The so-called discrete-time random- 
walk problem on a graph can be described by a matrix, 
therefore the distribution of its eigenvalue spectrum is 
much of interest. In this paper, we have studied the trap- 
ping problem on a class of regular treelike fractals with 
a deep trap located at the central node. The self-similar 



structure of these fractals can be characterized analytically 
in terms of their eigenvalues. We have derived a recur- 
sive relationship between eigenvalues of a relevant matrix 
for random walks on the fractals at different iterations by 
making use of renormalization, for which the computa- 
tional demands are very low. On the basis of the obtained 
recursion relation, we have characterized the complete set 
of eigenvalues and further determined their degeneracies. 
Moreover, we have provided a recursion relation between 
the smallest nonzero eigenvalues of the fractals at two suc- 
cessive iterations and obtained an approximate expression 
for the smallest nonzero eigenvalue. Finally, we have eval- 
uated the MFPT to the trap, showing that it is approxi- 
mately equal to the inverse of the smallest eigenvalue. 
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Proof of the equivalence of matrices Qg and Rg 

. — In order to prove that Qg = Rg, it suffices to show 
that their corresponding entries, denoted by {Qg)i.j and 
{Rg)ij, are equal to each other. Since both Qg and Rg 
are related to Pg+i (or Pg), we first determine the entry 
(Pg+i)jj of Pg+i before evaluating {Qg)i,j and {Rg)i,j. 

By the construction of the fractals, the element {Pg+i)i,j 
of the block matrix (see Eq. ^) can be easily found. Ac- 
cording to the properties of fractals, we have {Pa,a)i,j = 



5i^j, where 5i^j is the delta function defined as J^j = 1 if 
i — j, and Sij = otherwise. For {Pa,p)i,j and {Pp,a)i,j^ 
we distinguish two cases: if the nodes i and j are directly 
connected by an edge in Pg+i, then {Pa,i3)i.j = iPi3.a)i.j = 
— 1/ki, where ki is the degree of node i in Pg+i; otherwise 
{Pa.i3)i.j ~ {Pp.a)i,j ~ 0. Finally, for the diagonal block 
matrix P^,/3, each entry is given in Eq. ([3|). 

Notice that Qg is relevant to the inverse of the matrix 
Ag+i — P/3,0- So wc should determine (Ag+i — Pi3.p)~^. 
RecaU that P^_^ is block diagonal, then (Ag+i — Pp^p)^^ is 
also block diagonal with each block being (Ag+i — B)^^. 
Since our goal is to calculate {Qg)i,j, as will be shown, it 
suffices to find the entry of (Ag+i ~ B)~^ at its first line 
and first column, which is denoted by 0. By definition, 
we have 

(Ag+i - B)* 



(Ag+i - B)- 



dct(Ag 



B) 



(31) 



where (Ag+i — B)* is the complex adjugate matrix of 
(Ag+i — B). It is easy to verify that 



det(A 



9+1 



-P) = (Ag+l-l) 



m+1 



™(Ag + l - 1)^ 



which yields 

(A^+i - 1)^ 



e = 



+ 2) (A 



9+1 



1) 



dct( Ag+i -B) 2 + (m + 2) Ag+i (Ag+i - 2) ' 



(32) 



(33) 
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We are now in a position to determine {Qg)ij and 
{Rg)ij. We distinguish three cases. 

Case 1: i = j. In this case, it is easy to show that 
{Rg)i,j = X and 

{Qg\, =l + k, (^-1) e (-;^) =X= {Rg\j. 

(34) 

Case II: i ^ j , and at the same time nodes i and j are 
not directly Unked to each other in Fg. In the case, it is 



obvious that {Qg)ij = {Rg)i 



0. 



Case III: j but nodes i and j are directly connected 
in Fg. In this case, {Rg)ij = — y(— 1/fci) = y/kt and 

Thus, we have proved that Qg ^ Rg. 
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